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Abstract 

Based on a data sample of 60 pb~^ collected at the ^p{3770) resonance with the CLEO-c detector 
at CESR, we present improved measurements of absolute branching fractions for exclusive 
semileptonic decays into K''e'^u, tt'c^u and K*~e'^u; and the first observation and absolute 
branching fraction measurement of — > p~e~^u. 



* Submitted to the 32"^^ Int crnational Conference on High Energy Physics, Aug 2004, Beijing 
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I. INTRODUCTION 



The quark mixing parameters in the Cabibbo-Kobayashi-Maskawa (CKM) matrix are 
fundamental constants that must be determined from experiments. Semileptonic D meson 
decays are of great physics interest because they are relatively simple to handle theoretically. 
The decay matrix for semileptonic D meson decay decouples into a weak current (describing 
the Wlv vertex), and a strong current (for the Wcq vertex) that is parameterized through 
form-factor functions of the invariant mass (g^) of the W exchanged. The form factors 
can not be easily computed in quantum chromodynamics (QCD) since they are affected 
by significant nonperturbative contributions. That is the main source of uncertainty in the 
extraction of the CKM matrix element from the simple decay processes. Precise experimental 
measurements are needed to guide theoretical progress in this area. Charm semileptonic 
decays allow measurements of the form factors and CKM matrix elements Vcs and Vcd- 
Using Heavy Quark Effective Theory (HQET) or lattice gauge techniques, the measured 
charm form factors can be related to those needed to interpret h u transition and the 
measurement of the CKM matrix element Vub- 

While measurements of form factors in all charm exclusive semileptonic decays are im- 
portant, those in pseudoscalar-to-pseudoscalar transitions are the easiest to perform. The 
differential decay rate for the exclusive semileptonic decay D Piu {P stands for a pseu- 
doscalar meson) with the electron mass effects neglected can be expressed as yy]: 



dV 



where Gp is the Fermi coupling constant, is the four-momentum transfer squared between 
the parent D meson and the final state meson, is the momentum of the pseudoscalar meson 
in the D rest frame, and Vcq' is the relevant CKM matrix element, either Vcs or Vd- f+{<f) 
is the form factor that measures the probability that the flavor changed quark q' and the 
spectator quark q in Fig. n]will form a meson in the final state. The corresponding branching 
fractions can be obtained from 

r dV 

B{D ^ P£i^) =rDy< J dq^— (2) 



II. DETECTOR AND DATA SET 

The data sample used in this analysis was produced at the '?/'(3770) resonance with the 
Cornell Electron Storage Ring (CESR) and collected with the general purpose CLEO-c 
detector The total integrated luminosity of this data set is 60 pb~^ 

The CLEO-c tracking system covers about 95% of the 47r solid angle with a large cylindri- 
cal drift chamber and an all stereo wire vertex chamber. These tracking devices also provide 
specific-ionization (dE/dx) measurements for particle identification. The electromagnetic 
calorimeter is made of 7800 cesium iodide crystals with a total solid-angle coverage of 95% 
of 477. The calorimeter is crucial for electron identification and provides excellent efficiency 
and energy resolution for photons, yielding a typical mass resolution for 7r° reconstruction 
around 6 MeV/c^ (standard deviation). The Ring Imaging Cherenkov (RICH) detector 
covers 80% of the An solid angle and provides superb particle identification capability. 



3 



FIG. 1: Feynman diagram for charm meson semileptonic decays. 



III. EVENT RECONSTRUCTION AND SELECTION 

For this analysis we take fuU advantage of the unique kinematics of DD production 
at the '?/^(3770) resonance, just above charm threshold. We first select events in which a 
meson is fully reconstructed in one of the nine hadronic final states: i^~7r+, K^-k^ti^ , 
K~7T'^TT^TT^ , K~n~^7r~^7T~ , KsTT^, Ksi^^T^~ ■, KsTT^T^'Ti^, '7r+7r~7r°, and K~K^ . Charge con- 
jugation is implied here and throughout this paper. Within these tagged events we select 
the subset in which the meson decayed semileptonically to a specific final state. The 
efficiency-corrected ratio of the event yields gives the absolute branching fraction for the 
exclusive semileptonic decay mode. This branching fraction is independent of the luminosity 
of the experiment and benefits from the cancellation of many systematic uncertainties. This 
selection procedure provides an exceptionally clean sample of semileptonic decays that is 
ideal for determination of form-factor parameters and CKM elements, although our current 
data sample is not yet adequate for these determinations. 

Full details of the hadronic event reconstruction and tag selection can be found in a 
separate paper The selection is based on two variables. /S.E is the difference between 
the beam energy and the energy of the fully reconstructed candidate. M/j is the beam- 
constrained mass of the candidate, which is defined as Mn = \J — |p_Dp; where E^, is 
the beam energy and Pd is the measured momentum of the candidate. Fits to the beam- 
constrained mass distributions for candidates are shown in Fig. |21 Multiple combinations 
have been eliminated by selecting the candidate with the minimum value of |Ai?|. The signal 
component in these fits consists of a Gaussian and a bifurcated Gaussian, which is found to 
adequately describe initial-state radiation (ISR). The background component in these fits 
is presented by an Argus function Js^. The tag yields and efficiencies for selecting tags are 
given in Table IH 

The fully reconstructed meson serves as a tag and results in great suppression of 
backgrounds. We then identify an electron and a set of hadrons recoiling against the tag, 
to reconstruct the semileptonic decay side of the D^D^ system. Electron candidates are 
required to have | cos^| < 0.9, where 9 is the angle between the electron direction and the 
beam axis, with momenta greater than 200 MeV, the minimum to reach the Csl calorimeter 
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FIG. 2: Fits to the beam-constrained masses for different fully reconstructed decay modes. 
The signal is described by a Gaussian and a bifurcated Gaussian to account for the initial state 
radiation. The background is described by an Argus function j5||]. 



TABLE I: The tag yields, efficiencies, and signal yields and efficiencies for D — > K e'^v. 



Tags: -D° tag yields etag signal yields esignai 



K-7r+ 10183±112 

K'TT+TT^ 17208±160 

K-TT+TT^TT^ 4475±208 

K-TT+TT+TT- 14515±162 

KsTT+TT- 4629±101 

KsTT+TT-TT^ 4403±240 

KsTT^ 1502±43 

7r-7r+7r° 2533±176 

K-K+ 928±35 
all tags 



;65.35±0.12)% 245.7±15.3 

33.53±0.07)% 407.7±20.4 

16.10±0.06)% 106.6±10.6 

46.93±0.09)% 310.3±18.0 

37.05±0.15)% 100.0±10.2 

18.51±0.09)% 119.2±10.6 

27.92±0.22)% 26.6±5.3 

39.91±0.19)% 63.3±8.0 

;56.08±0.38)% 22.4±4.9 



(42.69±0.21)% 3.69±0.23 
(22.07±0.18)% 3.60±0.18 
(10.86±0.13)% 3.53±0.39 
(29.80±0.20)% 3.37±0.20 
(24.25±0.19)% 3.30±0.34 
(12.83±0.19)% 3.91±0.41 
(19.19±0.17)% 2.58±0.52 
(26.67±0.19)% 3.74±0.54 
(37.04±0.21)% 3.66±0.81 



1405.1±38.5 weighted ave. 3.52±0.10 
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for energy measurement. Electron candidates are selected with criteria that rely mostly on 
the ratio of the energy deposited in the Csl calorimeter to the measured momentum (E/p), 
dE/dx and RICH information. Extensive studies of efficiencies and misidentification rates 
have been performed using CLEO-c data for this electron identification procedure. Charged 
pions and kaons are identified using dE/dx and RICH information. For charged tracks with 
momenta less than 0.6 GeV, consistent dE/dx information (within 3 standard deviation 
(cr)) is required for pion or kaon candidates. For charged tracks with momenta greater 
than 0.6 GeV, additional requirement on RICH, if available, is imposed for pion or kaon 
candidates. K~ and tc^ candidates are combined to form K*~ candidates in K*~e^v 
decays. We require the invariant masses of the K*~ candidates to be within 100 MeV of the 
K*~ nominal mass In case of multiple K*" candidates, we select the candidate whose 
mass is closest to the K*~ nominal mass, vr" and ti^ candidates are combined to form p~ 
candidates in —>■ p^e^v decays. We require the invariant masses of the p~ candidates to 
be within 150 MeV of the p~ nominal mass In case of multiple p~ candidates, we select 
the candidate whose mass is closest to the p~ nominal mass. We reconstruct exclusive 
semileptonic decays into K~e^v, Tr'e^u, K*~e^v and p~e^v. No additional charged track 
is allowed to exist, besides those of the tagged and the semileptonic side of the other 
D^. The excellent particle identification of the CLEO-c detector helps in this reconstruction. 
The unique kinematics of threshold production provide additional and very powerful means 
to reject background from misidentified and missing particles. The difference between the 
missing energy and missing momentum in a event, U = E^iss — Pmiss, would peak at zero if 
the event is correctly reconstructed due to the undetected neutrino. 

In Fig. 121 we present kaon momentum versus U for signal — » K^e^v and the back- 
ground from D*^ TT'e^u reconstructed as — > K~e^v without hadron identification 
using MC simulation. It clear that the —>■ K~e^v signal and —>■ 7r~e^z/ background 
are well-separated even without kaon identification, thanks to the kinematics of threshold 
production. The separation decreases as the kaon momentum increases. 

Similarly, in Fig. |3J we present pion momentum versus U for signal —>■ Ti~e^v and the 
background from —>■ K~e^v reconstructed as — >■ T[~e^v without hadron identification. 
Again, the signal and the dominant background are well-separated. Therefore, tight hadron 
identification is not necessary for the semileptonic decay (signal) side. 



IV. RESULTS 

The absolute branching fractions of semileptonic decay is given by 

g -^signal /^signal 

-^tag/f tag 

-^signal is the number of the D^D^ events with one fully reconstructed and the other 
reconstructed through its corresponding exclusive semileptonic decay mode, egignai is the 
efficiency for constructing the fully reconstructed (tag) and the exclusive semileptonic 
decay of the Similarly, A^tag is the number of the fully reconstructed (tag) events 
observed, and etag is the efficiency for fully reconstructing the . Therefore, esignai/^tag is 
the efficiency for reconstructing the exclusive semileptonic decay in the presence of a 
found tag. 
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FIG. 3: Kaon momentum vs. U = -Emiss — Pmiss from D — > K~e'^v signal MC. The background 
from — > 'n~e'^v reconstructed as — > K~e'^v is overlaid and normalized by ^j^j^o^x-^e+u) ~ 
from PDG 0]. The right plot shows the projection onto the U axis in linear and log scales. 




FIG. 4: Pion momentum vs. U = Emiss—Pmiss from 'K~e'^u signal MC. The background from 

— > K~e'^u reconstructed as i:~e'^v is overlaid and normalized by Q^j^tiZllc-e+u) ~ 

from PDG [H]. The right plot shows the projection onto the U axis in linear and log scales. 
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In Fig.|3l we present the U distributions from data and MC for the selected D° K~e^v 
events with different tag modes. The comparison shows good agreement between the data 
and MC. Fig. IHl shows the decomposition of the U distributions from MC for all tag modes 
combined. The fit to the data is shown in Fig. ^| The K~e^v yields from various 

tag modes are given in Table EJ The efficiencies and branching fractions are also listed 
in Table HI for different tag modes. 
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FIG. 5: Comparison oiU = -Emiss — Pmiss for D K e^v between the data and MC. 



In Fig. [71 we present the U distribution from data and MC for the selected — » tt e+z/ 
events with all tag modes combined due to limited statistics. The comparison shows rea- 
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FIG. 6: Decomposition oiU = -E^miss — Pmiss for K~e~^u from MC for all tag modes combined. 

sonable agreement between data and MC. There is a clear excess in the — >■ 7r~e"'"z/ signal 
region (near zero) in data which is in good agreement with MC prediction. The background 
decomposition is also shown in Fig. [7| We fit the U distribution with a Gaussian, which 
corresponds to the vr^e+z/ signal near zero, another Gaussian which corresponds to 

the —>■ K^e^v background on the right side, and a 2nd order polynomial function. The 
2nd order polynomial is determined from MC simulation. The fit to the data is shown in 
Fig. El We observe 109.1±10.9 events. The signal efficiencies and branching fractions are 
given in Table ITTl 

In Fig. |S[ we present the comparison between the data and MC for K*-e+iy. The 

comparison shows good agreement between data and MC. The background decomposition is 
also shown in Fig. |H1 We fit the U distribution with a Gaussian and a 2nd order polynomial 
determined from MC simulation. The fit to the data is shown in Fig. E| The yields and 
efficiencies are given in Table ITTl 

The comparison of the U distribution for the data and MC for p~e^v is shown in 

Fig. El The background decomposition is shown in Fig. El We fit the U distribution with a 
Gaussian which corresponds to the p^e^v signal, another Gaussian for backgrounds 

mainly from K*^e^i>, and a 2nd order polynomial function determined from MC 

simulation. The fit to the data is shown in Fig. ^1 The yields and efficiencies are given in 
Table HD 

To determine Vcs and Vcd, we need to study the g^-dependent differential decay rate. 
Energy-momentum conservation is used to determine in the laboratory frame 

E\Y = Ef, — Ep, 

Pw = -PD - Pp, 

= E'^-\p^\'. (4) 

Here Ew and pw are the energy and momentum vector of the iu system or the virtual W. po 
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FIG. 7: U = Smiss — Pmiss for the selected — > vr e'^u events (left), and the background compo- 
nents (right). 



TABLE II: Efficiencies and yields for 
tag modes. 



vr 



^e'^u, — > K* e^v and —fp e^v using different 



Signal Modes 
Tags 



K*- 



e ' V 



p e' u 



K-TT+TT+TT- 
KsTT+TT- 

KsTT^ 

vr-vr+vrO 



(48.06±0.30)% 
(24.93±0.26)% 
(12.47±0.20)% 
(33.28±0.29)% 
(25.77±0.27)% 
(14.18±0.15)% 
(21.45±0.25)% 
(28.98±0.27)% 
(40.14±0.29)% 



(14.34±0.21)% 
(6.97±0.15)% 
(3.51±0.11)% 
(9.61±0.17)% 
(7.99±0.16)% 
(3.95±0.11)% 
(6.28±0.14)% 
(8.83±0.17)% 
(12.33±0.19)% 



(18.13±0.22)% 
(8.90±0.16)% 
(4.52±0.12)% 
(11.62±0.19)% 
(9.45±0.17)% 
(4.94±0.13)% 
(7.39±0.15)% 
(10.64±0.18)% 
ri5.14±0.21)% 



Yields 


109.1±10.9 


88.0±9.7 


30.1±5.8 


S(xl0~3) 


(2.46±0.25) 


20.69±2.28 


1.89 ±0.36 



is the momentum vector of the tagging D meson. Ep and pp are the energy and momentum 
vector of the hadronic system in the other D semileptonic decays. 

In Fig. ^2 we present distributions for K~e^u, — > ir^e^u and — K*~e^v 

from data. These distributions are from data directly, without any efficiency correction. 
Our q^ resolution is about 0.025 GeV^ or smaller which is over a factor of 10 better than 
CLEO III which achieved a resolution of 0.4 GeV^ This huge improvement in resolution is 
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FIG. 8: U distribution for the selected K* e^u events (left) and the components of back- 

grounds for — > K*~e~^h' (right). 




FIG. 9: U distribution for the selected p e^v events (left) and the decomposition of back- 

grounds (right). 

due to the unique kinematics at the ?/'(3770) resonance. The distribution for p^e^v 
from data is not shown due to the very limited statistics. Detailed study and measurements 
of form factor and CKM matrix elements will be performed with more CLEO-c data. 



11 



500 1—^ ~ ~ ~ 1 ~ ~ ~ ~ 1 ~ ~ ~ r 




>— ^^Tt-^-T -^rr-^-^— i^-^^ -i— J q ^^^^-^,-.;, , , 

-0.25 -0.125 0.00 0.125 0.25 -0.10 0.00 0.10 0.20 0.30 



U = E,;33-P„;33 (GsV) U = E^i33-p„.33 (GGV) 

FIG. 10: Fits to U = Smiss — Pmiss distributions for K~e'^u, 7r~e+i/, K*~e~^h' and p~e"*"z/, 

with the other fully reconstructed. 

V. SYSTEMATIC UNCERTAINTIES 

We have considered the following sources of systematic errors which may affect our results: 
charged track and vr'' finding, the veto on extra tracks, background normalization and shape, 
MC statistics, incomplete modeling of radiative corrections, model-dependent form factors, 
hadron identification, and electron identification. The contributions of these systematic 
errors are discussed below. 

From studies of CLEO-c data and MC, the systematic error for track finding efficiency is 
estimated to be 3% per charged track, 4.4% per 7r° and 3% per Kg. The systematic errors 
in particle identification are 2% for electron and 1% for pion or kaon. MC statistics varies 
from 1% to 3%, depending on the tag modes. To test the potential effect of uncertainties in 
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FIG. 11: Uncorrected distributions for —^K e^u, — > vr e'^v, ^ K*- 



radiative corrections, we have generated signal MC with and without ISR, and have found 
that the efficiencies vary within 1%. Using the selected K~e^v events, we find 0.51% 

events have extra tracks in our data, 0.41% events in our MC. Therefore, we assign a 0.5% 
systematics uncertainty due to the veto on extra tracks. To estimate the systematics from the 
backgrounds, we vary the background functions by ila, taking the correlation of parameters 
into account. We find the background normalization and shape contributes a systematic 
error of 1.1%, 3.1%, 2.9% and 5.3% for R-e+i^, n-e+u, K^'e^v and 

p^e^u, respectively. Systematic uncertainties on the reconstruction efficiencies from 
our current knowledge of form factors are small (< 2%) due to the uniform acceptance of 
the CLEO-c detector over large solid angle and momentum range. 

Most of these systematic errors are estimated from initial study of the first CLEO-c data. 
We expect these systematic errors to be reduced in the future, with more detailed study and 
more CLEO-c data. 

In Table ITTTl we summarize the contributions from systematic errors discussed above. The 
systematic errors from MC statistics is weighted with tag yields. We add these sources in 
quadrature to obtain the total systematic error for each decay mode. 



VI. SUMMARY 



In summary, we present improved measurements of exclusive semileptonic decays using 
the first CLEO-c data. The absolute branching fractions are listed in Table HVl 

Our results for — > K^e'^u and K*~e^v are consistent with those from the 

PDG [6]. Our B{D^ vr^e+z/) result is lower than that from the PDG 0. The ratio 
BiDoZK-e+% is closer to the CLEO III result of (8.2 ± 0.6 ± 0.5)%, while lower than the 
PDG value 0. The improved precision in these branching fractions is consistent with the 
expected performance of CLEO-c with initial data sample of 60 pb~^ at the ^/'(3770). The 
full integrated luminosity of ~3 fb^^ can be expected to provide great improvement in our 
knowledge of a much larger set of and exclusive semileptonic branching fractions as 
well as decay form factors and the CKM matrix elements Vcs and Vd |3 . 

We gratefully acknowledge the effort of the CESR staff in providing us with excellent 
luminosity and running conditions. This work was supported by the National Science Foun- 
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TABLE III: Systematics for four semileptonic decays (in %). 
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Li fXCrvlllg 


u 


u 


u 


U 


fn Ti 1 n (T 
/I iillClllig 








A. A. 


EID 


2 


2 


2 


2 


FiJJ 


i 


-1 
i 


-1 
1 


-1 
i 


extra track 


0.5 


0.5 


0.5 


0.5 


MC statistics 


<1 


1.1 


2.2 


1.9 


backgrounds 


1.1 


3.1 


2.9 


5.3 


ISR 


1 


1 


1 


1 


Form factors 


< 2 


< 2 


< 2 


5.0 


Yields 


1 


1.9 


1 


2.7 


Total 


7.0 


7.8 


8.9 


11.2 



TABLE IV: Absolute branching fraction measurements of the exclusive semileptonic decays, in 
comparison with PDG 6]. The uncertaities are statistical and systematic, respectively. 



Decays B_ PDG 

^ K-e^v (3.52±0.10 ± 0.25)% (3.58 ± 0. 
L>o ^ vr^e+i/ (0.25±0.03 ± 0.02)% (0.36 ± O.uoj/o 
Y)Q _^ K*-e^v (2.07±0.23 ± 0.18)% (2.15 ± 0.35)% 
L>o p~e+v (0.19±0.04 ± 0.02)% none 
gg;— (7.0 ±0.7 ±0.3)% (10.1 ±1. 
g%^'!r/r'*"-e+1) (9-2 ±2.0 ±0.8)% 



dation, the U.S. Department of Energy, the Research Corporation, and the Texas Advanced 
Research Program. 



[1] B. Grinstein, N. Isgur and M. B. Wise, Phys. Rev. Lett. 56 (1986) 258; F. J. Gilman and R. 
L. Singleton, Phys. Rev. D 41 (1990) 142; K. Hagiwara, A. D. Martin and M. F. Wade, Nucl. 
Phys. B 327 (1989) 569. 

[2] CLEO Collaboration, Y. Kubota et ai, Nucl. lustrum. Methods A 320, 66 (1992); T. Hih, 
Nucl. lustrum. Methods A 418, 32 (1998); D. Peterson et al, Nucl. lustrum. Methods A 478, 
142 (2002); M. Artuso et al, Nucl. lustrum. Methods A 502, 91 (2003). 

[3] CLEO-c Collaboration, CLNS-01/1742. 

[4] CLEO-c Collaboration, CLEO CONE 04-10, ICHEP04 ABSll-0775. 

[5] ARGUS Collaboration, H. Albrecht et al, Phys. Lett. B 241 (1990) 278. 

[6] Particle Data Group, K. Hagiwara et al, Phys. Rev. D 66, 010001 (2002). 



14 



[7] CLEO Collaboration, G.S. Huang et al, CLNS 04/1876, CLEO 04-06, CLEO CONF 04-14 and 
ICHEP04 ABSll-0780. Submitted to Phys. Rev. Lett. 



15 



